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Abstract 

We introduce a reduced model for a real sector of complexified Ashtekar 
gravity that does not correspond to a subset of Einstein's gravity but for which 
the programme of canonical quantization can be carried out completely, both, 
via the reduced phase space approach or along the lines of the algebraic quan- 
tization programme. 

This model stands in a certain correspondence to the frequently treated cylin- 
drically symmetric waves. 

In contrast to other models that have been looked at up to now in terms of 
the new variables the reduced phase space is infinite dimensional while the 
scalar constraint is genuinely bilinear in the momenta. 

The infinite number of Dirac observables can be expressed in compact and 
explicit form in terms of the original phase space variables. 
They turn out, as expected, to be non-local and form naturally a set of count- 
able cardinality. 

1 Introduction 

Ashtekar's variables ([I]]) simplify the algebraic structure of the constraint equations 
of general relativity so tremendously that one can solve various problems of classical 
and quantum gravity that were simply unfeasible before in terms of the old ADM 
variables. In particular, a number of model systems could be solved to more extent 
or even completely in these new variables. 

The model systems which could be solved completely up to now in terms of the new 
variables (e.g. 0) to our knowledge lack from not having in common at least one 
of the following features of general relativity : 

1) the reduced phase space is infinite dimensional (i.e. the number of Dirac observ- 
ables) and 
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2) the scalar constraint is bilinear in the momenta. 

Here we introduce a new model which is an algebraic plus Killing reduction of 
complexified Ashtekar gravity. In order to account for the fact that Ashtekar's con- 
nection obeys a non-trivial reality condition, we also impose reality conditions that 
bring us 'as close as possible' to the cylindrically symmetric waves (0) which form 
a genuine subset of Einstein gravity whereas our model defines then a real sector of 
complexified gravity which is manifestly no subset of Einstein gravity. 
However, our model has both features 1) and 2) and we are thus able to attack 
one of the major problems of general relativity, namely to find a complete set of 
(preferredly canonically conjugate) Dirac observables for a theory with an infinite 
number of degrees of freedom and non-trivial dynamics. 
The results of this paper are as follows : 

In section 2 we define the present model and show that it is not quite possible 
to formulate the cylindrically symmetric waves (which means to project our model 
of complex gravity to the real sector corresponding to Einstein gravity) in terms 
of Ashtekar's variables after we carried out the reduction. We then impose reality 
conditions that bring us at least 'as close as possible' to the cylindrically symmetric 
case. In particular, the first- and the diagonal part of the second fundamental form 
are real while the off-diagonal part of the latter is purely imaginary. 

In section 3 we carry out the symplectic reduction (|4[]) of our model. 
Quite surprisingly, it is straightforward to solve the typically very complicated scalar 
constraint before solving the diffeomorphism constraint. As usual for Hamiltonians 
or scalar constraints bilinear in the momenta, there are several 'sectors' of the con- 
straint surface. One of them corresponds to degenerate metrics, the other to non- 
degenerte metrics (the latter defines the 'physical' sector of our model). 
The most complicated step is then to find a complete but not overcomplete set of 
diffeomorphism invariant observables. This turns out to be feasible and the resulting 
functionals of the basic phase space variables can be given in explicit and compact 
form (the observables given in H are rather complicated and untractable objects 
due to the fact that the mode system of Bessel functions is involved). 
As expected for a diffeomorphism invariant field theory, the Dirac observables come 
out to be naturally smeared and nonlocal. In contrast to non-diffeomorphisminvariant 
field theories one does not have to employ a complete system of mode functions on 
the hypersurface to write the Dirac observables in this form. The reduced phase 
space thus has naturally a countable cardinality. We argue that this is always true 
for a spatially diffeomorphism invariant field theory. 

Section 4 is devoted to quantum theory. 
From the reduced phase space point of view this is already trivial since the phase 
space derived in section 3 can be cast into a form so that it has the structure of a 
cotangent space over a real, infinite dimensional vector space. 

The application of the algebraic quantization programme forces us to make use of 
functional integral techniques. 

The two resulting Hilbert spaces turn out to be unitarily equivalent. 
As to be expected for a field theory with constraints quadratic in the momenta, 
the final set of Dirac observables mix configuration and momentum variables with 
respect to the original polarization (here the Ashtekar polarization). 



2 



While it is important to derive the Dirac observables, they are typically hard to 
interprete. Here the framework of deparametrization (||) helps to define 'evolving' 
Dirac observables which allow for a direct interpretation in terms of metric variables 
and to leave the limits set by the 'frozen' formalism. 

We extend this formalism to field theories and apply it to our model. We thus 
arrive at an interpretation of the theory although the meaning of the constant Dirac 
observables remains veiled. 

Finally, in section 5, we discuss the loop representation (|7|]) of our model. 
First we define a set of loop variables that stand in bijection with the phase space 
variables underlying the connection representation. This implies that we can ex- 
plicitely express constraints, symplectic form and Dirac observables in terms of loop 
variables. These expressions are quite complicated compared to those obtained for 
the connection representation and that allows for two possible conclusions : 
Either it seems to be doubtful whether the loop variables are tailored for the prob- 
lem to construct the reduced phase space or 

Killing reduced models are unsuited to simulate the situation in the full theory as 
far as the loop representation is concerned. 

We summarize what has been learnt in dealing with the present model, valid for 
the full theory of Einstein gravity. 

2 Definition of the model 

We assume the reader to be familiar with the Ashtekar formulation of canonical 
gravity (0). In the sequel a,b,c,.. will denote tensor indices and i,j,k,.. SO(3) 
indices and we apply the abstract index formalism. 

We impose the following three conditions on the phase space variables of the com- 
plexified Ashtekar formulation of gravity : 

1) Killing- reduction : all fields are Lie-annihilated by a set of two commuting, hy- 
persuface orthogonal Killing fields. 

Since the Killing fields commute they are surface forming. Since they are hypersur- 
face orthogonal these surfaces lie orthogonal to a foliating vector field. Hence we 
find three global independent coordinates x,y and z such that the Killing fields are 
given by d x , d y . Then all fields depend on z (and the time variable t) only. 

2) algebraic reduction : We set to zero the following components of the basic fields 
of the Ashtekar formulation 

A\ = A\ = Al = A\ = E{ = E* = ES = E y 3 : = . (2.1) 

Up to now we are on the phase space of the so-called Gowdy models ([§). We will 
assume, for simplicity, the spatial topology of the initial data hypersurface to be 
that of the three-dimensional torus in order to avoid technical difficulties associated 
with boundary conditions. 

One more condition brings us close to the axisymmetric spacetimes. 

3) algebraic reduction : The Killing fields are orthogonal. 
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This last condition forces the densitized triad to be of the form 



(Et) = 



I E x cos(P) -EP sin(/3) 

E x sm(P) EPcosiP) j (2.2) 
V E z 



and in order to have as many configuration variables as momentum variables we 
require in analogy to (2.2) the Ashtekar connection to reduce to 

A x cos(a) —A y sm(a) 
(A\) = | A x sm(a) A y cos(a) j (2.3) 
A z 

where all variables are generally complex up to now. 

As far as the triads are concerned, provided we restrict them to be purely real, we 
easily see that we have already a real diagonal 3-metric (compare the appendix). 
Let us see whether we can also implement the reality condition 

Re(A i a - rj = 0, T\ = spin-connection (2.4) 

which leads to the Einstein phase space. In order to do that we first have to compute 
the spin-connection. The calculation is deferred to an appendix. The result is given 
by 

/ -r x sin(/3) r y cos(/3) \ 
(P )= T x cos(P) T y sm(P) (2.5) 

V o or,y 

where 

1 E y E z , 1 E X E Z , , 

Tx = 2E~y (y ~E x ~" > ' Ty = ~2E x ~ (y ~E~y~" > ' Tz = ~^ ' (y2 ' 6 " > 

The problem is already obvious : both, T l a and A l a have five nonvanishing com- 
ponents, giving rise to five reality conditions induced by (2.4), but these compo- 
nents are coordinatized respectively by only four independent coordinates (namely 
A x , A y , A z ,a and E x , E y , E z , (3). The explicit form of these reality conditions is 



given by the requirement that the following quantities {iK l a := iK ab E\ f/y det(g) = 
A\ — Y\ where K ab is the extrinsic curvature) 

%K X := A x cos(o;) + T x sin(/5), iK x := A x sin(o;) — T x cos(/3), 

iKy := — AySm(a) — T y cos({3), iK y := A y cos(a) — T y sm((3), 

iK z := A z - T z (2.7) 

are imaginary. Equivalently 

%K X E\ = E y (A x sm(a - 0) - T x ), iK l y E x = E x (-A y sm(a - (3) - F v ), 

%K X E X = E X A X cos(a - (3), iK l y Ef = E y A y cos(« - (3), 

iK\E z = E Z (A Z - T z ) (2.8) 

are imaginary. 

We will now show that these are five functionally independent equations whence our 
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reduction is incompatible with (a subset of) the real form of complexified gravity 
that corresponds to the Einstein phase space : 

We will immediately show that the Gauss constraint of this model is given by 



g = (E z y 



sm(a - [3)(A X E X + A y E y ) . 



Let 5 := a — (3. Then, using the Gauss constraint we derive that 



tan(<5) 



(2.9) 



E X [A X cos(<5)] + Ev[A y cos(5)] 



which is weakly imaginary since the square brackets in the denominator are accord- 
ing to (2.8). Assuming that the imaginary part of 5 is finite, this implies unam- 
biguously that 5 itself is imaginary : 5 — i£ where £ is real. Then cos(5) = ch(£) is 
real whence again from (2.8) we derive that A x , A y are imaginary. Finally we have 
7 := A z + a' = (A z — T z ) + 5' = (A z — T z ) + i£' which displays 7 and A z — T z as 
(weakly) imaginary quantities. 

We still have to satisfy the two first conditions in (2.8). First of all, again using the 
Gauss constraint, we find that 



where we have used the identity (E z )' — T X E X + T y E y = (which is due to the 
covariant Constance of the triads) and the reality of the triads. This demonstrates 
that the real part of the first quantity in the first line of (2.10) vanishes if and only 
if the last line in (2.10) is zero, at least on the Gauss-constraint surface. However, 
since A x , A y , 5 are imaginary, the imaginary part of both expressions is already zero 
so that at least one of the fields A x ,A y ,5 would be expressible completely in terms 
of the triads and therefore the symplectic structure would become degenerate. This 
furnishes the proof. 

Since we are mainly interested in playing with a model that captures some of the 
features of general relativity rather than describing the cylindrically symmetric case,[] 
we take the following viewpoint : 

We want to preserve at least the reality of the diagonal part of the extrinsic curvature 
(the extrinsic curvature is diagonal for the cylindrically symmetric waves). Since 
the metric is manifestly diagonal, we can impose the simple reality condition that 



are all purely imaginary. Now, looking at (2.8), we see that the extrinsic curvature 
will be symmetric upon imposition of the Gauss constraint but it will develop an 

lr The cylindrically symmetric case is not of physical importance anyway because an axisymmet- 
ric matter distribution is quite improbable 




Re(- A X E X sin(6) + E X T X ) 



(2.10) 



a. — f3, A z — T z , A x , A. 



y 



(2.11) 
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imaginary off-diagonal part in the course of evolution. The only case when it is real 
is when the off-diagonal part vanishes and we are then back on a genuine subset 
of the phase space of the cylindrically symmetric model (this subset is a fix point 
under evolution). The diagonal part is real. The metric is diagonal and real by 
construction so we managed to impose reality conditions that are 'as close' to the 
cylindrically symmetric possible. 

These reality conditions are also preserved under evolution as we will show shortly 
and that raises another problem : Due to the evolution law of full general relativity 

q ab = 2NK ab + D {a N b) (2.12) 

an initially real, diagonal metric becomes complex and nondiagonal when inserting 
our K ab above. This is apparently a contradiction ! However, there is a mistake in 
this argument : the processes of deriving the field equations and reducing down the 
degrees of freedom cannot be expected to commute if the phase space of the reduced 
model does not lie entirely in that of the full theory as it is the case for our model. 
One might ask whether it is possible to recover the cylindrically symmetric case 
by making use of all SO (3) degrees of freedom (the three Euler angles rather than 
only one, a and (3 for A l a and Ef respectively). However, if that was true one could 
always go to the gauge (2.2), (2.3) and since the reality structure of Einstein gravity 
is SO(3) covariant (since SO(3) is a real Lie group) one would obtain a contradiction. 
Thus, it is not possible to treat the cylindrically symmetric case while making use 
of a corresponding reduction on the Ashtekar phase space. 

Nevertheless, this model has all features in common with full general relativity 
(except for the modified reality conditions) so that it serves as an interesting new 
testing ground for the quantization programme of full quantum gravity as we will 
now show (actually this is also the viewpoint that one takes in dealing with 2+1 
gravity as a toy model for 3+1 gravity : the connection of 2+1 gravity is purely 
real) . 

3 Complete symplectic reduction of the model 

It will turn out that the symplectic reduction programme can be carried out com- 
pletely, that is, with respect to all constraints. 

After integrating over the (finite) range of the coordinates x and y, the reduced 
action becomes we neglect a trivial prefactor coming from this integration) 

S= f dt f dz[-iAlE? - (iA l Qi - iN a V a + n\c\ (3.1) 

where it is understood that we have to substitute for A l a , Ef the reduced expres- 
sions (2.2) and (2.3). Here Qi, V a , C are, respectively, the Gauss-, Vector- and 
scalar constraint and A 4 , N a , N are the corresponding Lagrange multipliers, called 

respectively the Gauss scalar potential, shift vector and N := ^Jdet(q)N the lapse 
function. 

By plugging the formulae (2.2) and (2.3) into the expressions for the constraint 
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functions we obtain 



Qi := V a Ef:=d a Ef + €i jk AiE^ 

= S i3 [(E z Y - sin(a - (3){A X E X + A y E y )} =: 5 l3 G 

= S z [cos(a - mA x )'E x + (A y )'E y ) - A Z (E Z )' - a' sm(a - P)((A X )'E X + (A y )'E y ] 

= S z a [cos(a - P)((A X )'E X + (A y )'E y ) - 7 (£ 2 )' + olQ\ =: 8 z a V 

C 1 ^ Z7* Z? a J?b 

G .- e ijk l< ab E j h k 

= 2[A x A y E x E y + {{A X )'E X + (A y )'E y )E z sin(a - /3) 

+7(A r £ :E + A y E y )E z cos(a - =: 2C (3.2) 

where is the curvature of the Ashtekar connection, is the Levi-Civita symbol 
and we have abbreviated 7 := A z + a'. Further we have absorbed the trivial factor 
of 2 in the last line into the Lagrange multiplier of the scalar constraint. 
Now we insert the reduced form of our basic coordinates into the symplectic potential 
and obtain up to a total differential (note that spatial integrations by parts do not 
contribute boundary terms because the torus has no boundary) 

iS [d t ] = dz[A z E z + A X E X cos(a - (3) 

+A y E y cos(a - (3) - a(A x E x + A y E y ) sin(a - (3)} 
= dz[(A z + a)E z + A X E X cos(a - (3) + A y E y cos(a - (3) + 6cQ\ 

^[ 7 n 7 + Ajr + A y u y + «n a ] . (3.3) 



IT, 

Note that all configuration variables on the Gauss-reduced phase space are imaginary 
while the momenta are all real. This is a very nice reality structure. 
We next write the constraints in the so defined canonical coordinates. For the Gauss 
and vector constraint this is easy : 

Q = n Q , 

V = (A X )'U X + (A y )'Tl y + a'Tl a - 7 (rF)' (3.4) 

while for the scalar constraint we have to do some more work. 
We have directly from the definitions of our canonical coordinates 

(IP)' - IT* 

tan(« -0) = 4rr rW" (3-5) 

v 1 J A x n x + A y W v ; 

so that we have managed to express the combination a — (3 in terms of canonical 
variables. 

Next, we use the trigonometrical identities 

. 2 tan 2 2 1 



sin = 77, cos 



1 + tan 2 ' 1+tan 2 
and have unambigously 
IPIP 

C = A x A y — — - + {{A X )'W + {A y )'W)Wt^{a-f3)+ 1 {A x W + A y W)W 

cos z (a — (3) 
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A x A v U x U y (l + [ n " I 2 ) + ((A X )'U X + (A V )'W)W 

(wy - ip 

V / I _Y A TTX i /I TT3/MT7 



a x ii x + A,n> 



+ 7(AJP + 4TF)ir 



1 -) 2 A x A y u x uy([(wy - irf + [a^if + ^irf ) 



aji x + 4,nw' 

+((A x yu x + (A y yuy)W[(wy - u a \[A x u x + A y u y \ + 7(4,11* + ^n^n 7 } 
= : W + ^n^ ' (3 ' 6) 

We now rescale all configuration variables except for a by a factor of i so that 
these become real too (e.g. — iA x will be called A x again. Then the reduced action 
becomes 

S = [ dt f dz[A x E x + A y E y + 7IT 7 + idAT - [iAU a + N X V - NC}} (3.7) 

with manifestly real constraint functions Q, V and C and Lagrange multipliers. This 
furnishes the proof alluded to in the previous section that the reality structure of 
the model is preserved under evolution. 

We will assume that the prefactor 1/(A X U X + A y Il y ) 2 never diverges or vanishes so 
that we can absorb it into the Lagrange multiplier and thus obtain a fourth order 
scalar constraint for this model. 

It is not trivial to check whether the constraint algebra closes. Before doing that, 
we will pass to the following equivalent set of constraints, obtained by subtracting 
a term proportional to the Gauss constraint from the scalar and vector constraint : 

G = IT\ 

v = (A x yu x + (A y yu y - 7 (ir)' 

C = A x A y WW(-((Wy) 2 + (A x W + A y W) 2 )+ 1 (A x W + A y WYW 

-i(A x yu x + (A y yuy)w(wy[A x u x + A y u y ] (3.8) 

and something amazing happens when adding the term II 7 (II 7 ) '[Axil 1 + Aj / Il y ]l / 
proportional to the vector constraint to the scalar constraint 

c = A x A y n x u y (-((wy) 2 + (A x n x + A y u.y) 2 ) 

- 7 n 7 [(n 7 ) / ] 2 [^n :r + A y w\ + 7(4,11* + A y wfw 

= [-((IF)') 2 + (A X U X + A y U y ) 2 ][A x A y U x U y + 7 n 7 (A c IP + A y U y )\ (3.9) 

i.e. the scalar constraint factorizes into two second order constraints with respect 
to the momenta ! They are given by 

d := A x A y U x U y + 7A c iriT 7 + ^A y Tl y TV and (3.10) 
C 2 := -((n 7 )') 2 + (A X U X + A y U y ) 2 . (3.11) 

By studying the transformations that these constraints generate, we obtain that 

A x , A y , 7, n x , IF, II 7 , Il a are 0(2)-invariant while a is an 0(2)-angle and that A x , A y , a, IV 

are scalars while E x , E y , IT*, 7 are densities of weight one, i.e. Q generates 0(2) 
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transformations while V generates diffeomorphisms on the torus. By the way, 7 
can be interpreted as the abelian version of a connection plus the Maurer-Cartan 
form &mc — +(dS)S~ 1 (which in one dimension for S = exp(a) boils down to 
Omc — ot'dx) and therefore is an 0(2)-invariant quantity. 

The only nontrivial Poisson-bracket to check is then again that between two scalar 
constraints (all constraints are 0(2) invariant and Diff(X) covariant). We have 

{C[M],C[N]}= Jjx fjyM{x)N{y){C x {x)C 2 {x),C x {y)C 2 {y)} 

= Jjx ^^M^iV^^^^^I^W^^^j + C!^)^^!^^),^^)} 
+C 2 (x)C 1 (y){C 1 (x), C 2 (y)} + C^C^C^x), C x {y)}\ 

= \jjx ljy(M(x)N(y)-M(y)N(x))[C 1 (x)C 1 (y){C 2 (x),C 2 (y)} 
+2d (x) C 2 {y){C 2 {x),C 1 {y)} + C 2 (x)C 2 {y){C x (x) ,C 1 (y)}} 

= Jjx Jjy{M{x)N(y) - M(y)N(x))C 1 (x)C 2 (y) 
[-5 >x (x,y)]2(Wy(x)W(y)(A x U x + A y Uy)(y) 

= J dx{M'N - MN^dC^YTfy^W + A y W) 

= C[{M'N - MN')({U r ) 2 )'{A x U x + A y W)} (3.12) 

i.e. the constraints form a first class algebra, open in the BRST-sense, and we can 
apply the framework of symplectic reduction (0|). 
We begin with the Gauss-constraint : 

We just have to pull back the symplectic potential to the IP = surface, the vector 
and scalar constraint as they stand in (3.8) and (3.9) are already reduced. 
As experience with other model systems showsQ, it turns out to be much more con- 
venient to first reduce the scalar constraint and then the vector constraint because 
the scalar constraint is not diffeomorphism invariant but only diffeomorphism co- 
variant. 

We have to distinguish between two possible sectors : sector I is defined by the 
vanishing of Ci, sector II by the vanishing of C 2 . 

We should mention that sector II corresponds to metrics that are degenerate in the 
z-direction : For our model the metric tensor takes the form 

EVE* E Z E X E X E\ 
q ab = diag(-^-, -^-) . (3.13) 

From the Gauss-constraint and the definition of U x , U y we infer 

{A X E X + A y E v ) 2 = -{(E*y? + [A X W + A y Wf and g = (3.14) 
so that we can solve E x ,E y , E z in terms of reduced coordinates 

(E x , E y , E z ) = (± - - W A - Jc 2 , ± - - W A - JC2, n 7 ) (3.15) 
v ' ; v A X H X + AyUv v ' A x Il x + AyW v ' ; y ' 

2 specifically, the author was looking at spherically symmetric gravity plus matter (|9j) 
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so that 

n^n 7 ifif ifif 

q ab = diag(± — , ±— , ± [AxTLx + AyTLv]2w C2) ■ (3.16) 

Note that the quotient W/W remains finite. 
This shows that sector II is the 'unphysical' one. 

3.1 Sector I 

We pass to new canonical pairs, 

(q x := HA X ),P X := A x W;q y := ln(A,),^ := A y W; q := ln( 7 ),p := 7 IF) (3.17) 

because then the scalar constraint Ci adopts the simple form 

pV +p y p + pp x = (3.18) 

which does not even make any ordering problems any more if one would pass to 
the quantum theory directly, that is using the Dirac approach. Moreover, it defines 
trivially an abelian constraint because it consists only of momenta. 
Note that we could have changed the polarization for the pair 7, II 7 in order that all 
configuration variables are scalars without destroying the simple form of the scalar 
constraint, but then we would leave the Ashtekar-polarization. Let us therefore 
proceed this way although this will equip the canonical coordinates with a somewhat 
awkward tensor valence. 

It is clear that one will now diagonalze this constraint because then one can solve 
the Hamilton- Jacobi equation by seperation of variables. We have 

C = \[{p x +P y + 2pf - (p x -p y ) - (2p) 2 ] (3.19) 

so that we define still another set of canonical pairs 

P°:=^(p x +py + 2p), pi:=±( p *- p y), P 2 := p (3.20) 

and conversely 

f = p° + p 1 _ 2P 2 , p y = P°-P 1 - 2P 2 , p := P 2 (3.21) 

which allows for the determination of the new canonical configuration variables via 
inserting (3.21) into the symplectic potential 

e[d t ] = J^dx[q x (P + P 1 -2P 2 )+q y (P°-P 1 -2P 2 ) + qP 2 ] 

= jf dx[P°j t (q x + q y ) + P 1 ^ - q y ) + P 2 |(g - 2q x - 2q y )} (3.22) 
from which we read off 

Qo = q x + q y , Qi = Qx- q y , Q 2 = q ~ 2q x - 2q y . (3.23) 
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Then the scalar constraint just says that (P°, P 1 , P 2 ) is a null- vector in a 3-dimensional 
Minkowski space at every point of E. There are two branches of general solutions 
to the associated Hamilton- Jacobi equation, 

S±\p»,P v -M = J s dz[±y/(p«)> + (p«) 2 Q + P u Qi + P V Q 2 ] (3.24) 

where, following the second reference of ||, p u = P ,p v = P 2 are the new momenta 
and 

q v := £ = ± , PU Q° + Q l 

5p u ^( p «)2 + (pvy 

q v := £ = ± . ^ Q° + Q 2 (3.25) 

are the new invariant (under the motions generated by the scalar constraint) con- 
figuration variables. 

One can get them also just by pulling back the symplectic structure to the scalar 
constraint surface via its embedding i into the phase space (compare ||). Inserting 
the solutions P° = ±J (P 1 ) 2 + (P 2 ) 2 into the symplectic potential yields (as before 



we do not display total differentials) 

{i*e)[d t }= [ \±Q^(p u ) 2 + ( Pv ) 2 + Q lP u + Q 2 p v ] 

= -/ [p"(± pU Q + Ql ) + f( ± P° Q + Q 2 ) 

= I \p^(± ^ =Q + Q 1 )+p v ^(± P 1^ =Q + Q 2 ()3.26) 

It remains to reduce the vector constraint. Resubstituting all the symplectomor- 
phisms we have carried out up to now, we obtain (note that for any two functions 
a and b holds 

b VL {a 2 + b 2 )\, {a 2 + b 2 Y , 2 



Va 2 + b 2J Wa 2 + b 2 ' 2Va 2 + b 2 2(Va 2 + 6 2 ) 3 

.IF. 



v = (A x yu x + (A y yw-^(pp 



= (q x )'p x + (q y )'p y -l( — y 

7 

= (qx)'p x + (q y )'p y + q'p - p' 

= (Qo)'p° + (Qi)'p 1 + (g 2 ) / p 2 - (p 2 )' 

= ±(g )V(p u ) 2 + (p v ) 2 + (Q^'p u + (Q*Yp v - (p v )' 

= (q u )'p u + (q v )'p v - (p v )' 

= (q u )'p u + (Hq w ))'q w p w - (q w p w )' 

= (qu)'p u - q w ip w )' (3.27) 
where we have carried out one more canonical transformation 

q v =: \n(q w ), p v =: q w p w . (3.28) 
Equation (3.28) displays q u ,p w as scalars and q w ,p u as densities of weight one. 
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3.1.1 Reality structure of sector I 

From their definition (3.20) it is clear that p x ,p y ,p as products of two real quantities 
are all real, whence also p u ,p v are real. 
From (3.17) we have that 

Qo = m(AA), Qi = M^r), Qi = Htj\ti) (3-29) 

and so we see that the reality structure of q u , q v depends on the sign of the arguments 
of the logarithms above. However, they are always real possibly up to an imaginary 
constant ±in/2 which drops out of the symplectic potential anyway if we assume 
the number of sign changes of the configuration variables for any instant of time to 
be of measure zero with respect to the measure dx on the hypersurface (hence, the 
set {x G £; lim a _>o(sgn(j4 x + s) — sgn(Aj — s) ^ 0} is such that d/dt\n(sgn(A x )) 
etc. are integrable where sgn(:r) := x/\x\ is the sign of x). Hence we will assume in 
the sequel that is real whence this is also true for q u , q v . This then implies that 
q w = exp(g„) is always positive. 

3.2 Sector II 

Now the scalar constraint C 2 can even be written as one of two branches of a 
constraint linear in momentum 

C± := A X U X + A y W ± (TP)' . (3.30) 

Quite similar as for sector I we define new canonical pairs 

(q x := ln(A x ),p x := A x U x ;q y := ln(A,),^ := A y Il y ; q := 7 ,p := IF) (3.31) 

and obtain 

C ± :=p x +p y ±p f . (3.32) 
The general solution of the Hamilton- Jacobi equation is obviously given by 

S ± = j^dx[{p v )'q x + {p w )'q y =F (p v +p w )q] = Jjx\p"{ T q ~ (&)') + P w (Tq ~ (q y )')} 

(3.33) 

whence we read off for the new momenta 

(p v Y=p c , ( P w y = P y (3.34) 

and for the new configuration variables 

q v ■= Tq - (qx)', q w ■= Tq - (q y )' (3.35) 

which are both densities of weight one so that the reduced (with respect to the scalar 
constraint) vector constraint becomes 

v = -q v (p v y - q w (p w y = (qu)'p u - q w (p w y ^-m 
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where we have defined still another canonical pair by 



Qv = 



: (q u )', q u := -(p")' . 



(3.37) 



That this defines indeed a symplectomorphism is obvious from the following short 
calculation (spatial surface integrals vanish due to our choice of the spatial topology) 



3.2.1 Reality structure of sector II 

Again we have that p x ,p y ,p are all real which then is true via (3.34) also for p u ,p v . 
With the same assumption as for sector I we obtain that q v , q w are real and finally 
via (3.37) also for q u ,q w . This time, however, q w is not necessarily positive. 



3.3 Reduction of the vector constraint 

We come now to the main problem for a spatially diffeomorphism invariant field 
theory, namely to find a complete but minimal set of Dirac observables which are 
to be promoted to the basic set of quantum operators later. 

It is always easy to find (classical) expressions that Poisson-commute (even strongly) 
with the diffeomorphism constraint - the recipe is the following : 

1) From the transformation law of the basic field variables under diffeomorphisms 
derive their tensor nature. 

2) Construct scalar densities of weight one. 

3) Integrate them over the initial data hypersurface. 

These objects are then already invariant under diffeomorphisms that are connected 
to the identity (in case of an asymptotically flat topology in general only if the 
diffeomorphism tends to the identity rapidly enough at spatial infinity). This is 
obvious from the transformation property of a scalar density s of weight one under 
an infinitesimal diffeomorphism (generated by a vector field £ a ) 8s = d a (^ a s). 
However, the scalar constraint is not at all so easy to solve since it can be viewed as 
a dynamical constraint. Hence, it is most important to keep the algebraic structure 
of the scalar constraint as simple as possible. 

Now, since the scalar constraint is a density of weight two, it has a nontrivial trans- 
formation law under diffeomorphisms. Accordingly, solving the vector constraint 
before solving the scalar constraint will most likely complicate rather than simplify 
the algebraic structure of the scalar constraint because it will become a non-local 
object and, moreover, cannot be reduced because it cannot be written in terms of 
diffeomorphism invariant quantities. 

The conclusion of all this is that the scalar constraint is to be solved before the dif- 
feomorphism constraint for a field theory with a scalar constraint which is at least 
quadratic in the momenta. This is what we were able to do up to now. 

Now we are going to derive a sequence of results which generalize in an obvious 
manner (although this will become combinatorically much more difficult) to higher 




(3.38) 
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dimensions. 



Lemma 3.1 The set of all integrated scalar densities constructed from the configu- 
ration variables provides for a (possibly overcomplete) system of coordinates of the 
diffeomorphism reduced configuration space. 



1) Any local object constructed from the configuration variables has a nontriv- 
ial transformation law with respect to infinitesimal diffeomorphisms. Accordingly, 
diffeomorphism-invariant objects are non-local, i.e. they can be expressed as inte- 
grals over (parts of) S of local objects. 

2) Were these integrals not over all of £ then they were not invariant because the 
diffeomorphisms have their support everywhere except for the boundary of E. 

3) Only scalar densities of weight one are invariant when integrated over S. 

4) One could consider objects of the form 



i.e. K(x\, ..,x m ) is some integral kernel. However, any such integral kernel can be 
approximated arbitrarily well (in the given norm) by linear combinations of tensor 
products fi(xi)..f m (x m ) (for each k the different functions fk will belong to a base). 
Thus, the integral of K can be written as a sum of products of single integrals which 
displays K as a derived object. 
This proves the lemma. 
□ 

Let {q l }iei denote this set of Dirac observables where the labelling is with respect 
to an index set I whose cardinality is to be obtained. Let us further assume that all 
these observables are algebraically independent. Any linear combination 



of these observables is again an observable where the pi are, in general, complex 
numbers (the reality conditions on the pi match those on the q l so that there are as 
many independent momenta as configuration variables). Then, provided that the 
set of the q l is complete in the sense that every diffeomorphism invariant object 
constructed from the configuration variables can be expressed by them, the above 
functional S is the Hamilton- Jacobi functional relative to the diffeomorphism con- 
straint and the pi play the role of the integration momenta (recall the formalism in 
the second paper of 

We will now define an overcomplete set of diffeomorphism invariant functionals of 
the configuration variables for our model and then systematically obtain a minimal 
subset. When stripping off the technical methods that are particular for this model, 
one will get an insight for what is necessary for the full theory. 
According to lemma (3.1) what we have to do is first is to construct all possible 
scalar densities s that can be built from q u , q w . We first need two technical lemmas 



Proof : 




(3.39) 



(3.40) 
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Lemma 3.2 1) In one dimension, tensors of valence (r,s) are scalar densities of 
weight s-r. 

2) In one dimension, a scalar density of weight one s gives rise to a s-compatible 
connection T according to 

r := pn(S)]' . (3.41) 

Proof : 
1) 

This lemma is most easily proved by comparing the corresponding transformation 
properties of tensors T a - 6 .. under diffeomorphisms x a — > x a = x a — £ a 

5T a - = C.I-'- = CI'"- - V • ,,..C, ~ ■■ + '/" • r..C,, + •• (3-42) 

which in one dimension reduces to 

ST = £T' + (s - r)i'T . (3.43) 

On the other hand, a scalar density s of weight n transforms as 

6s = Z a s >a + nst% (3.44) 

which again in one dimension reduces to 

§§ = £§' + n§£' (3.45) 

and that furnishes the proof of the first part of the lemma. 

Accordingly, in one dimension tensor valences are just characterized by their density 

weight. 

2) 

One simply builds the unique torsion-free connection compatible with a given met- 
ric by restricting the Christoffel formula T^ c = l/2q ad (q dbc + q dcb — qb c ,d) to one 
dimension. A metric in one dimension is, according to 1), just an arbitrary density 
of weight 2 corresponding to a symmetric tensor of valence (0,2) (i.e. a metric q). 
Since the inverse of a metric in one dimension is just one over this metric we have 
for the affine connection in one dimension T = l/2(q + q — q)/q = l/2[ln(g)]' or, 
by defining the density s of weight one through the equation q = s 2 , we have finally 
a connection defined by an arbitrary density s of weight one given by 

T := pn(5)]' . (3.46) 

Again, by just restricting the covariant derivative 

D c T a - b _. = d c T a - b _. + T a cd T d - b _. + .. - T d cb T a - d .. - .. (3.47) 

to one dimension we find that for a density T of weigt n we have 

DT — T' — nTT (3.48) 

and one can explicitely check that DT is a density of weight n+1, i.e. a density 
again with one more covariant valence. 

In particular Ds = s' — Ts = which shows that the connection is s-compatible. 
□ 

We are now ready to prove one of the main theorems of this section. 
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Theorem 3.1 An (overcomplete) system of scalar densities of weight one is given 
by 

q w f({(qu) (n) }neM) (3.49) 

where f is an arbitrary function of the infinite number of arguments indicated, M 
is the set of non-negative integers and for any scalar function s, denotes the 
scalar d n (s) where d is the operator d = l/q w d which is obviously anti- self- adjoint 
with respect to the measure dji = q w dz. 

Proof : 

It is clear that any scalar density s of weight one can be written as q w s where s is 
a scalar, simply define s := s/q w . The proof thus boils down to showing that the 
above function f in (3.49) is the most general scalar that we can build. 
Such a scalar will be the most general scalar constructed from q w ,q u and all their 
partial derivatives up to arbitrarily high order. This can be done in a covariant way 
only by contraction of tensors, that is in one dimension, by multiplying densities 
such that their overall weight is zero, according to lemma (3.2), first part. Hence we 
need to construct a covariant derivative such that the derived tensors, i.e. densities, 
transform as tensors, i.e. densities again. This can be done by constructing an affine 
connection in one dimension and here part 2) of lemma (3.2) comes into play. 
We are now able to take arbitrarily high covariant derivatives of q w , q u , for example 
with respect to the connection derived from s := q w . But then Dq w = because D 
is q w compatible which shows that the covariant derivatives of q w do not contribute 
in our construction and furthermore D n <f) = q™(l/q w D) n <f) = q™d n <j) where we have 
introduced the operator d = 1/ q w d. 

It remains to show that there is no loss in generality in choosing s := q w . 
Assume that we had chosen any other density s instead of q w in the construction 
of T then the covariant derivatives D of q w would not vanish, however then Dq w = 
sd(q w /s) = sD(q w /s) = —q w /sDs and s is again built only from our configuration 
variables and D (it must be a tensor and it cannot be constructed by using D because 
D is defined through s). Hence Dq w can be recast in the form of a function f of the 
type given in (3.49). 
This furnishes the proof. 
□ 

What we have achieved up to now is to obtain a sufficient number of observables q % 
which are functional of the configuration variables only. 

We now ask for possible redundancies among them. Two sources of redundancy are 
obvious : 

1) by doing integrations by parts it is possible to relate some of the q % and 

2) there might be topological identities which could arise, for example, via index 
theorems and thus affect the range of q % (for example in 2 dimensions the integrated 
densitized curvature scalar is an object of the form we are considering here, but 
the Gauss-Bonnet theorem states that it depends only on the genus of the compact 
hypersurface without boundary). 

Under the assumption that we can restrict to analytical f's in (3.49) it is possible 
to get rid of the first redundancy. What we are presupposing is that there exists an 
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expansion of f of the type 



f{{a m }meAf) = f{n k } II a k" ( 3 - 50 ) 

{n k }} k=0 

where := d k q u , the sum is over all sequences of (non- negative) occupation numbers 
rik and the coefficients of the monomials 

oo 

M {nk} ■■= n < I 3 ' 51 ) 

A;=0 

are real numbers because q u , q w are real as shown above. Hence, in our case the index 
set is given by the set of sequences, I = {{nk}} 3 i, and the (still overcomplete) set 
of q v s is given by the integrals 



[ dzq w Mi (3.52) 



whereas the /j are the conjugate momenta. 

We will devide out the first, combinatorical redundancy by bringing the q % into a 
standard form thereby introducing a terminology that one is used to from Fock- 
space or statistical mechanics. : 

First we observe that by doing an integration by parts the occupation numbers nk 
of the quantum states or energy levels k for particles of type q u may change but that 
the total particle number (of each particle type seperately if there are more than 2 
fields on the configuration space which it is not the case in our situation) 



as well as the total energy 



N:=J2n k (3.53) 

fc=0 



E:=J2 kn k (3.54) 

fc=0 



are preserved under such a process, simply because the number of factors of a (ir- 
respective of how often they are derived) and the number of derivations does not 
change. In fact, one could write the operation of doing an integration by parts in 
terms of annihilators and creators. 



Lemma 3.3 Consider now the class of (q l ) 's with the same total particle number 
N and total energy E. Then an algebraically independent genuine subset of q % 's in 
this class is given by those with indices i such that mo = l..N,mi = 0, the rest of 
the occupation numbers being arbitrary (subject to the condition that they lie in the 
class labelled by (N,E). 

Proof : 

In the first step we show that every member in the given class can be written as a 
linear combination of these g*'s and in a second step we show that these q l cannot 
be transformed into each other. 
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Consider the class (N, E) = (M, k) and we may exclude the trivial case M — 0. 
Then q % is of the form 

q l = J ' dfj. a kl ..a E -(ki+..+k M - 1 ) (3-55) 

where we do not care about occupation numbers, i.e. some of the k r may be equal. 
Let, without loss of generality, k\ be the lowest of the k r . If k\ = we are done, 
otherwise do an integration by parts. After that, q % is a linear combination (with 
integer coefficients) of q^s of the required form since all these g J have m = 1. 
By doing further integrations by parts if necessary for these g- 7 (since /c(a ) fc ~ 1 ai = 
d((a ) k ) can be integrated by parts) one can satisfy mi = 0,m > 0. This finishes 
the first step. 

Let now one of these 'reduced' q l, s be given. Then the set of its occupation numbers 
is characterizing, i.e. by doing an integration by parts, starting from m\ = 0, one 
necessarily picks up terms that have mi ^ from the differentiation of a since we 
already have m > 0. This furnishes the proof. 
□ 

From now on we will understand the index set I to be reduced, that is, we consider 
the 'Mandelstam-identities' between the original q l, s originating from integrations 
by parts to be taken care of by lemma 3.3 . 

As far as the second source of redundancy is concerned, we do not have any indication 
that any of the q l should have a non-continous range. 
There is still a more subtle source of redundancy : 

One could fix a gauge (a frame) and expand q u into a taylor series with respect to 
this frame where we have assumed these variables to be smooth functions on E. The 
freedom in choosing these functions is then captured in the choice of the values of 
their taylor coefficients. 

However, then one could just integrate out all of our reduced q % with respect to this 
frame and we would get complicated rational functions of these Taylor coefficients 
multiplied by real numbers coming from the integration of the powers of the frame- 
variable x over the hypersurface and depending only on the topology of E. 
Since any such rational function can be generated from the Taylor coefficients itself, 
the knowledge of the latter is sufficient and provides for the minimal subset of 
observables. 

The result of all this analysis summarized in the subsequent theorem. 

Theorem 3.2 Only the following q l 's are independent variables on the diffeomorphism- 
reduced configuration space (n is any non-negative integer) : 

q n ■= \ £ dfxq u d 2n q u = (-l) n ^ jf dfi[d n q u ] 2 . (3.56) 

Proof : 

There is a one to one mapping between these functionals and the Taylor coefficints 
of q u , possibly up to sign. For example, after fixing a gauge, e.g. q w (x) = g'{x) = 
1, d — d where g is a fixed scalar function, the q n allow to extract all the Taylor 
coefficients of q u (e.g. d n 4> starts with the n'th coefficient), possibly up to sign. 
□ 

Note that (3.56) are \M\ independent quantities, where \S\ denotes the cardinality 
of the set S, because they satisfy all the criteria of the previous theorem to be 
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algebraically independent and precisely for that reason we have b n / 2 = & n /2 = 
0. Accordingly, the number of degrees of freedom is naturally countable without 
employing a complete system of mode functions on E to express the true degrees of 
freedom in discrete form. 

Moreover, these observables have the nice feature that 

= q w d 2n q u (x), (3.57) 



Sq u (x) 



-\ 2 ^{-l) k {d k q u {x)){d^ qu {x)) . (3.58) 



8q w (x) 2 



k=i 



We will denote the momenta conjugate to q n by p n ). 
The associated Hamilton- Jacobi functional is 

oo 

S:=X>"<Z« (3.59) 

n=0 

and the old momenta are the functional derivatives of S with respect to the old 
coordinates. 



For sector I we have seen that q w is positive (and therefore serves as a genuine 
density for a measure on X in the the definition of q n ) whereas for sector II the 
range of q w is the whole real line. Since q u , q w are both real, (d n q u ) 2 is positive 
whence q n has half the or the full range of the real line for sector I or II respectively. 
p n on the other hand is generally real so that the pair (q n ,p n ) coordinatizes the 
cotangent space over the (positive) real line. 



Our, admittedly somewhat sketchy, construction generalizes in an obvious way to 
more dimensions, for every diffeomorphism invariant field theory. 



4 Quantum theory and determination of the phys- 
ical inner product 

4.1 Reduced phase space approach 

Quantum theory is straightforward. Keeping the polarization that we have ended 
up with, physical states depend only on the variables q n . 

It is important to see that the scalar constraint forced us to let the reduced configu- 
ration space (with respect to the scalar constraint) not only depend on the original 
configuration variables (before reducing) but also on the original momenta ! This 
mixing of original polarizations is typical for genuinely quadratic constraints with 
respect to the momenta and can be viewed as a source of difficulty in the Dirac 
approach to select physical states. 

The Dirac observables q n ,p n are both real where p n has infinite range and q n has in- 
finite or half infinite range for sector II or I respectively. Therefore, the unique inner 
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product that turns q n ,p n into self-adjoint operators is determined by the physical 
Hilbert spaces TC p h ys 



L 2 (X™ =0 R, A dq n ) or L 2 (x~ B+, f\ dq n /q n ) 



(4.1) 



ra=0 



n=0 



respectively (as follows from the group theoretical approach [10|). 



4.2 Dynamics of the model - deparametrization for an infi- 
nite number of degrees of freedom 

What we have obtained up to now is the complete set of Dirac observables which 
are time-independent. However, their interpretation is still missing. On the other 
hand, interpretation is usually straightforward for non-gauge invariant quantities. 
In order to arrive at an interpretation there are two possible strategies available : 
either one expresses an object that is 0(2) x Diff(E) invariant in terms of Dirac 
observables plus an additional gauge function t which labels the choice of foliation of 
spacetime and quantizes only the Dirac observables (||11||) or one uses the framework 
of time-dependent Dirac observables (deparametrization, ||). 
Let us try to apply the latter one. 

In the sequel we will only deal with the 'physically relevant' sector I. 
We are asked to cast the scalar constraint into a form that allows to apply an exten- 
sion to an infinite number of degrees of freedom of the framework of deparametriza- 
tion. As is proved in (||) the only change is that the time variable becomes 'many- 
fingered' (Tomonaga-Schwinger extension). 

In particular we are interested in the time-development of the quantities E x , E y , E z 
which determine the metric. 

The starting point will be equation (3.19). We take the quare root and restrict to 
the positive branch. Then the scalar constraint adopts the familiar form 

P°(z) - H(Q 1 , P 1 ; Q 2 , P 2 ){z) := P°(z) - ^[P\z)Y + [P*{z)\* (4.2) 

which is already deparametrized, although at each point x of the hypersurface seper- 
ately. In (4.2) we have defined an infinite number of 'Hamiltonians', H(x), that is, 
an 'energy density'. Note that, since H involves only momenta, there are no ordering 
problems. 

Following this framework, we construct the time-evolution operator 

U t := exp(-z / dz(Q - t)H) (4.3) 

and define any operator O (to be factor-ordered appropriately) on the set of physical 
states \I/ = Uij) (with respect to the scalar constraint) by 

(d t ^)[Q ,Q 1 ,Q 2 \ := uA=qMQ\Q 2 ] ■ (4-4) 

Now, in contrast to toy models for cosmology, physical states should also be dif- 
feomorphism invariant. We will therefore take the operator to be diffeomorphism- 
invariant, but not necessarily invariant under the scalar constraint. 
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We want to apply the above framework for example to the operator for the three- 
volume whose classical counterpart is given by 

V := J dz^j&et q = j dzVE x EvE z . (4.5) 

We will proceed as follows : We start from (4.2) and proceed to the Dirac quan- 
tization of the model. As for the scalar constraint, we have no operator-ordering 
problems and the complete solution of 

(-[P ] 2 + [P 1 ] 2 + [P 2 ] 2 mQ , Qx, Q 2 ] = (4.6) 

is formally given by a functional integral (positive energy solutions) 

tt[Qo,Qi,Q 2 ] = / [d 2 k](j)[k\ k 2 } exp(-i / dz[^{k l f + (k 2 ) 2 Q - k 1 Q 1 - k 2 Q 2 }) 
= expH JjzQoy/iP 1 ) 2 + (P 2 ) 2 )^[Q 1; Q 2 ] (4.7) 

where [d 2 k] is an appropriate limes of a cylinder measure in the usual sense (skele- 
tonized Lebesgue measure (fll2||) with a finite number of integration variables since 
£ is compact and the limes corresponds to the removement of the skeletonization) 
and (ptplk 1 , k 2 ] is a yet arbitrary functional. 

Since U t commutes with the vector constraint (ordered with the momenta to the 
right), diffeomorphism invariance is obeyed if and only if for an infinitesimal density 
£ of weight -1 holds 

1>[Qi + (QiYt, Q2 + (Q 2 )'£ + £] = 1>[Qi, Q2) ■ (4.8) 

Inserting, the integrand of the path integral becomes after doing an integration by 
parts in the exponential (up to first order in £) 

<f>[k\ k 2 } exp(* f dxiik 1 - (£k 1 )')Q 1 + (k 2 - (£k 2 y)Q 2 + i'k 2 }) . 

We change to new variables K l := k % — (£&*)', i = 1, 2 in the functional integral and 
obtain (up to first order in £) the condition 

j>[K x + K 2 + (K 2 0'] exp(* jf dzt'K 3 ) = cf>[K\ K 2 ] (4.9) 

where we have used the diffeomorphism invariance of the Lebesgue measure (the 
Jacobian is 1 to first order in £ : denote by h the spacing of the sceletonization, 
k = k 1 or k 2 , k n := k(nh), then the Jacobian is 

dk n d[K n - l/(2h)((£K) n+1 - (£K")n-i)] 



dK m dK m 

= 8 n . m — —(S nm -i — S nm+ i) =: J nm =: (1 + A{^)) nm (4-10) 
2n 

where the matrix A is trace-free and of first order in £, so det( J) = 1 + tr(A) + .. = 

1 + o(e)). 



21 



We now split ^[k 1 , k 2 } =: exp(i J s dzf(k 2 ))^[k 1 , k 2 ] and seek to determine the func- 
tion f according to the requirement that 

jf dx[f(K 2 + {UK 2 )') + £'^ 2 ] = jf ^/(iT 2 ) 

up to first order in £. We expand f and perform integrations by parts until only £ 
appears derivated with respect to z. Then a sufficient condition for the last equation 
to hold is that 

*(l + ^>-,-0. 
The solution of this ordinary first order differential equation is given by 

/(A; 2 ) = [c-ln(A; 2 )]A; 2 

where c is a real integration constant. 

This means that ip is diffeomorphism invariant when k 1 , k 2 transform as densities of 
weight one. Again we are therefore lead to the functionals 

#\* 2 ] = /({?"}), where q n = J^{<T{k 2 /k 1 )) 2 (4.11) 

with dfj, = dzk 1 , d = l/k l d. 

Diffeomorphism invariance up to first order implies diffeomorphism invariance with 

respect to the component of the identity of the diffeomorphism group. 

We are now in the position to define the operator V. 

By making use of the definitions we can express (4.5) in terms of Q^, P^ 

px/y 



EX/V = "I 7 -J(P x +P y Y- TO 2 hence 

(p° - p 2 ) 2 - (p 1 ) 2 



^/det(q) = exp(-Q°) 



2 (po _ p2) 

^4( P o _ P 2)2 _ [( ex p(-(Q 2 + 2Q Q ))P 2 )'} 2 (4.12) 

and we face severe factor ordering problems when promoting the integral over this 
function to quantum theory. 

The action of this operator on physical states is given by (we assume that we order all 
momenta to the right handside in the formal power expansions of the non-analytical 
terms in (4.12)) 

VMQo,Qi,Q2] = U t J[d 2 k}exp(iJ^dz[Q 1 k 1 + Q 2 k 2 + [c-\n(k 2 )]k 2 ])x 

k 2 ) 2 -{k 1 ) 2 

J^yjk 2 

(4.13) 



jf dx exp ( - 1) Vl _ , v / 4 ( A;2 ) 2 -[( ex P(-(^ + 2t))A; 2 )'] 2 



which, when acting with the evolution operator (which, by the way, is unitary for 
Qo — t =const. with respect to the inner product (4.1) : since the integrated Hamil- 
tonian commutes with both constraints, it is physical, that is, well-defined on TC p h ys , 
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and since it is classically real it promotes to a self-adjoint operator by a suitable 
choice of ordering) results in a horrible object which we cannot work out explicitely. 
The point of this subsection was to show how deparametrization works in principle 
for field theories and how one can arrive at an interpretation. It also shows that 
most of the evolving operators which one would intuitively think of as observables 
are probably not well-defined on the physical Hilbert space (here in the connection 
representation). 



5 Loop variables 

The concluding section of this paper is designated to the issue of using (traces of) 
Wilson-loop observables for the quantization of or even for the classical treatment of 
gravity. Because our model system is completely integrable one can ask and answer 
all kinds of questions that have been raised in the literature. 

In particular we ask if the operators that are usually defined for the loop represen- 
tation in the literature (@) are well defined at all, that is, if they make sense when 
applied to physical states. Of course, as they stand they will not leave the physical 
Hilbert space invariant and thus are not even expected to have a well-defined action 
on physical states. However, it will be interesting to see how the Dirac observables 
look in terms of loop-coordinates, in particular the constraints, what their Poisson 
structure is and how one can define a 'loop-transform' into the loop representation 
(0), from states that have been defined for the connection representation. 
Also, one could construct loop-operators as certain 'time-dependent' Dirac observ- 
ables. 

We begin by recalling that one defines the following 'basic' loop variables when deal- 
ing with the loop representation : 

The T° variables are defined to be traces of the holonomies (parallel transport op- 
erators) 

U a [A}(x) := V ex P ( [ A) (5.1) 

around closed curves (loops a) with base point x in the fundamental representation 
of SU(2) (the Ashtekar variables are complex- valued whence one they lie in the Lie 
algebra of complexified SU(2,C)=SL(2,C); although the gauge group of gravity is 
really SU(2), the loop variables are even invariant under SL(2,C)) 

T a [A] := \tx\UM\] (5-2) 

where V denotes path-ordering and A := A l a dx a T i , Tj a basis in the Lie algebra 
su(2). We will choose T; L := — i/2a iy the latter being the usual Pauli matrices. Note 
that the trace of the holonomy is independent of the base-point on the given loop. 
The ^-variables are defined to be 

T*[A,E](x) :=te[E a {x)U a [A\{x)\ • (5-3) 

As one can show, they form a closed Poisson-algebra and can thus serve as a starting 
point for the group-theoretical quantization scheme. Let us explore what happens 
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with these variables for our model. 

The (pull-back under a section of the principal SU(2)-bundle of the) connection 
reduces to 

A(x,y,z) = A x (cos(a)Ti + sin(a)T2)dx + A y (—sm(a)Ti + cos(a)T2)dy+A z T 3 dz (5.4) 

and the gauge transformations generated by the Gauss-constraint for our reduced 
configuration variables are those under the U(l) subgroup of SU(2) generated by t% 
(even its complexification) : 

- dSS- 1 + SAS- 1 where S = exp(Ar 3 ) (5.5) 

can be written as in (5.4) only that a, A z are replaced by a + A, A z — A'. Under a 
gauge transformation, the holonomy and the triad E = E^Tid a transform covariantly 

E{x) -> S(x)E{x)S-\x) and U a(x) -> S(x)U a (x)S-\x) . (5.6) 

whence the T°, T 1 variables live on the Gauss- reduced phase space. 
It is convenient to eliminate the angle a by doing an U(1,C) rotation with gauge 
parameter A = —a (note that a is complex so this transformation is really an 
element of U(1,C)). Then the connection and the triad adopt the reduced structure 

A(x,y,z) = A x (z)ndx + A y T 2 (z)dy + ^y(z)r 3 dz (5.7) 
E(x,y,z) = [U x (z)d x - Ey(z)sm(a- (3)d y ]r 1 + [W(z)d y - E x (z)sm(a- (3)d x ]r 2 
+W(z)r 3 d z . (5.8) 

Note that both quantities do not have any x,y dependence ! 

We will now start to find suitable loops to capture the full U(l) invariant informa- 
tion needed in order to actually coordinatize the Gauss-reduced phase space. Three 
loops are obvious : they are those lying entirely in one of the three ^-factors of the 
three-torus. We will call them a(x),a(y) and a(z) respectively. A short calculation 
yields for the corresponding T°,T l, s (the range of the three coordinates x,y and z 
is taken to be the interval [0, 1] with endpoints identified) : 

T a(x) (z) = cos(^A x (z)) (5.9) 
T a{y) {z) = cos{\A y {z)) (5.10) 



2 

T a{z) = cos( 



\f Q d Z1 {z)) (5.11) 

T% x) (z) = -Il*(z)sm(±A x (z)) (5.12) 

T: (z) (z) = -Uy(z)sm(^J o 1 dz 1 (z)) (5.14) 

T ^ - (^S^ W^ (5 - l5) 

T: (z) {z) = Tl (z) {z) = T: {x) {z) . = T* [y) {z) = (5.17) 
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Note that from the four quantities ^,^,(2), fi G {x,y} we can already ex- 
tract the two gauge invariant canonical pairs A^, IF locally while the two vari- 
ables T5 z -j (z) , To{z) allow to recover only II 7 locally whereas the integrated quantity 
J E dz r y(z) is a global one. This has two consequences : 

1) The information contained in the two 'off-diagonal' variables Tjf, v Jz), n 7^ v is 
already redundant and so we can neglect them, 

2) We need a further T° variable in order to get hand on 7(2). 

Before doing that we show that the variables defined above form a closed (clas- 
sical) *-Poisson algebra (the Loop-variables for the Ashtekar phase space are not 
closed under complex conjugation at least it is not obvious how to express T a = 
|tr(Pexp(/ Q [2r — A])) in terms of T°, T 1 where T is the spin-connection). 
We have proved in the sections before that A x ,A y ,^f are purely imaginary while 
IP, IF, IF are purely real. This implies the *-relations 



T a{a) (z) = T a(a) (z) and T% a) {z) = -T% a) {z) (5.18) 

i.e. the T°'s are purely real, the T^s purely imaginary, and also that the range of the 
T° are the real numbers larger than or equal to one while the T x 's have the whole 
imaginary axis as their range. This opens the chance for a cotangent topology of 
the associated phase space and therefore a group theoretical quantization scheme. 
Another consequence is that the variables T^qa are not quite sufficient to recover 
the connection up to gauge transformations : we have the degeneracy 

T aM [A] = T aM [-A) (5.19) 

which means that these loop variables are not seperating on the moduli space A/Q of 
SL(2,C) connections modulo gauge transformations for our model when restricting 
oursek'lves to the above defined class of loops. This is in analogy to the known 
degeneracies of these loop variables as coordinates for A/Qf\ 



Next we compute the Poisson-structure which turns out to be quite similar to the 
one obtained in the first reference of (H) when restricting ourselves to the z-direction 



{T <a) {z),T a{h) {z')} = (5.20) 

1 

1 0,1 r 1 , , , , 1 



{T ctM (z)X m (/)} = U^z^sm^A^z)) = ^M(V')(1 - [TvUz)] 2 ) 



(5.21) 

{T* M (z),T^ z) (z')} = ^sin 2 (^ i 7 (^)^) = i(l-[T (T{z) ] 2 ) (5.22) 
{T^(zlT: {u) (z')} = {2^(20^,0*0} = (5-23) 
{T* (z) (z),T* {z) (z')} = ~[T* {z) (z)-T* {z) (z')]T a(z) (5.24) 

and we could even have written the right handside of all the equations as linear 
combinations of single T°'s or , T^s with distribution- valued coefficients at the price 
of extending the so defined loop algebra to all such loops with any winding number 



3 This result implies, for example, that any toplogy on A/Q which is based on the use of the 
loop variables (traces of the holonomy) displays A/Q as a non-Hausdorff space (Q) ! 
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(with respect to an orientation on the three torus) by making use of the Mandelstam 
identities 

T a [A}T p [A] = i[T ao/3 + T aof3 -i) (5.25) 

(H) where o is the product in the loop group (with respect to a given base point). 
This Mandelstam identity becomes here just the trigonometrical identity 

cos(a) cos(/3) = ~[cos(a + 0) + cos(a — 0)\ (5.26) 

but we refrain from introducing the winding number because we will have to deal 
with higher order polynomials of loop variables anyway. 

Note incidentally that (5.21) implies that the the conjugate variable to T^,^ is given 
by 2arth(T .( At) ) ! 
Remark : 

It is clear that the Hamilton- Jacobi functional S := cT CT ( z ) solves all constraints for 
the degenerate sector for any constant c. Thus it turns out that the loop variables 
serve as 'good' coordinates only for (part of) the degenerate sector as for the full 



theory (Hl6| ). 



We now have to extract j(z). One way would be to make use of the so-called area- 
derivative (||15||), whose relevant components for the already defined loops become 
(the area-derivative in the (ab)-plane is denoted 5 a b) at the point with coordinate z 
(x,y are irrelevant) 

1 If 1 

SxyT a ( g ) = --A x (z)Ay(z) sm(- dzj(z)) (5.27) 

S yz T a{x) = ~A y (zh(z) sin(~ A x (z)) (5.28) 

S zx T a{y) = -^(z)A x (z) sin(± A y (z)) (5.29) 

but since the area-derivative is a certain limit of a difference between one of the 
above loop variables and an extended one (the area derivative in the (ab)-plane with 
respect to a given loop at one of its points is obtained by attaching an infinitesimal 
loop lying in the (ab)-plane to the given loop at the point of interest, computing 
the trace of the holonomy for this modified loop and the original one, taking the 
difference and dividing by the area element) we actually leave the above algebra and 
look at an infinite number of new loops. We would like to have only as many loop 
variables as reduced canonical variables. Accordingly, we are looking for a better 
suited object. 

After some trial, one of the most simple objects seems to be an eyeglass loop (j7|) 

obtained by drawing two loops with equal orientation at z\ and z 2 in x or y direction 

connected by a line in z-direction back and forth between z\ and z 2 . 

We will denote this loop by o{xz) or <r{yz). The computation of the associated T° 

reveals 

11 1 1 

Ta^izt, z 2 ) = cos(-A„(zi)) cos(-A x (z 2 ))-cos( / dzj(z)) sm(-A x (zi)) sin(-A x (z 2 )) 
2 Z J zi 2 2 

(5.30) 
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and thus allows to extract j(z) by taking the derivative with respect to, say, z 2 = z. 
Fix a value zq on the z-circle and write 



T a ( xz ){z) := T a{xz ){z ,z) . (5.31) 

We observe that 

T a(xz) (z + 1) = [T a{x )(z )} 2 - (2[T CT(Z) ] 2 - 1)(1 - [T a(x) (z )] 2 ) (5.32) 

whence 



1 — T a ( xz )(zo) 
\\2(l-[T a(x) (z )Y 



(5.33) 



so that the information contained in T a i z \ is also contained in T a t xz \ (z) (if one uses 
also T a (jj\(z)) (note that even its sign can be determined from (5.23) as the T c r a \(z) 
are all positive). Therefore we will consider T a ^ as redundant in the following al- 
though we will use it as an abbreviation for (5.33). Also the analogon T a i yz \ of T a ( xz \ 
is not needed. 

We now compute the extended loop algebra : the only non-vanishing brackets of 
our new element with the rest are (up to a sign) 

{T a{xz) (z),T* {z) (z')} = - X [z ,z]( z ') sin ( f Z -r(t)dt) sin(i f dtj(t)) 

J zq Z JO 



sin(-A x (z )) sin(-A x (z)) = -X[z , z ](z') x 



~2 2 1 J 1 ' 

y/(l - [T„ {z) mi - [T ct( ^ )] 2 )(1 - [T„ {x) {z)Y) - (T a(xz) (z)-T aix) (z )T aix) (z)y} 

(5.34) 



x - 

(5.34) 

1 

+5(z', Zi )[-T: (x) (z ) +T: {x) {z)T a{xz) {z)]) (5.35) 



{T a{xz) (z) } T* {x) (z>)} = -(5(z',z )[-T: (x) (z) +T: {x) {z )T a{xz) {z)\ 



Unfortunately, the right handside becomes quite complicated and it would seem to 
be better to work with an overcomplete set and apply the algebraic quantization pro- 
gramme ([fL7|1). In our case that would boil down to the problem to find loop- variables 
that depend only on the configuration variables, contain s\n{\A^{z)), sin(/ 1 dz^(z)/2) 
and yield brackets that have analytical right handsides for the corresponding Poisson- 
brackets, but that turns out not to be straightforward without regarding the whole 
set of loops. 

We now come to set up a bijection (up to signs) between the Loop-variables defined 
so far and the canonical variables determined in section 3. 

We have already defined the T's in terms of A^, IP; 7, II 7 . Here is the local inversion 



A^z) = 2axcos(T a ( lx - ) (z)) (5.36) 
T M (z) 

lP(z) = ; " M (5.37) 



7 (z) = [arcos( , r »(')fo) W*) ~ ^(**)(*) )y (5 . 38) 
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TP(z) = T j^ Z) (5.39) 

1 ~ 2(l-[T ff(l) (* )] 2 ) 

where the prime denotes derivation with respect to the z-argument only. The z- 
derivation considered as a difference of Loop- variables does not lead out of the class 
of variables that we have defined and so we are fine. 

(5.36)-(5.39) defines a (local) point transformation which is always a local symplec- 
tomorphism. Accordingly when plugging the above formulae into the symplectic 
potential, we are actually able to determine the momenta conjugate to the three 
variables T CT(ju) , T a{xz) ! 

After some tedious algebraic manipulations, an integration by parts, using that 
Is 1 l{z)dz written in terms of T°'s is a spatial constant and making use of the fact 
that quantities frozen at z are also spatial constants we derive the momenta con- 
jugate to 

T a ( x )(z), T a ( y )(z), T a ( xz )(z), the basic loop variables that we are using, respectively 
as 

" T a{x) (z)T a{xz) (z)} 



x 



y/(l - [T^)P)[(1 - [T a(x) (z )}i)(l - [T a{x) {z)Y) - (T a(xz) (z) - T aix) (z )T a{x) (z)y 

x 1 T Z(x)(z) , , 

i - (T a{x) (z)r i-[T„ {x) (z)r 1 • } 



i-\T a(y) {z)r 



(5.41) 



y/(l - [T a{z) n(l - [T a(x) (z W)(l - [T a(x) (z)Y) - (T a(xz) (z) - T a{x) (z )T a{x) (z)y] 

(5.42) 

and the momentum conjugate to T a ^(z ) is the integral over S 1 of expression (5.40) 
only that the roles of T a ( x )(z ) and T a ^(z) are interchanged. The reason for the 
appearence of this global degree of freedom is that we fixed the point z to be an 
observable by hand. 

Next we look at the constraints expressed in terms of loop-variables. From their 
definition it is obvious that T a ( p ) is a scalar, and since (5.38) shows that the density 
7 is a spatial derivative of a function of our three loop variables, we find that also 
T a ( x ,z) must be a scalar (alternatively, this follows from its definition). T^) (^o) is 
a constant and its conjugate, as an integral over a density, is also diffeomorphism 
invariant. 

Thus, necessarily the vector constraint must be 

(T aM )>S» + (T a{xz) yS (5.43) 

where we have denoted the corresponding conjugate momenta in (5.40), (5.41) by 
S^, S. Of course, it is also obvious from the definitions that and are 
densities and a scalar respectively. 

It turns out to be more convenient to make direct use of the inversion formulas 
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(5.36)-(5.39) to obtain the explicit expressions 



rpl rpx rpl rp x 

y _ or "( g ) _|_ a (y) a (y) 



+ [arcos(^== y — r _ n )] / ; ] (5.44) 



1 - (Ta(x)) 2 1 - (T ff(! ,) 

) rp r 

Sff-iiy 2 )' J n ^-*2u 

( , = 4 ^ arco < T -w) T ^) arcos ( T -fe)) T o'(j/) 



1-[T" ,] 2 



| r arcos ( T CT fa))^ fa ) arcos(T CT(:E) )T; (;c) ^ 



x [arcos( ^M Z )~ T <^ y ^kL=] (5.45) 

' [^)] 2 )(i-KV)] 2 ) 9 _7E^ 

2 1-[T°. ,]2 

where all expressions are to be evaluated at z except for those with a superscript 
which are to be evaluated at z . 

The result looks much more complicated than in the connection representation. 
From the point of view of the loop representation one would have to take equations 
(5.44) and (5.45) and apply it (after a choice of ordering) to a loop state ([0]) ac- 
cording to the Rovelli Smolin loop representation of T°, T 1 . In order to do that one 
would have to expand the trigonometrical formulas involved in both constraints and 
one would obtain an expression involving an infinite sum of loop states. This is a 
very discouraging result. 

It is clear how to obtain the observables of the theory : simply substitute the loop 
variables via (5.36)-(5.39) into the expressions for the observables given in section 
3. These formulas are too long and it is not worthwhile as to display them here. 



6 Conclusions 

Let us briefly summarize what we have learnt by studying this model system : 
First of all we worked out a number of ideas and methods to reduce a spatially 
diffeomorphism invariant field theory which might be useful to keep in mind for 
future work. In particular, our methods do not suffer from the problem of the 
existence of a bijection between a frame on the hypersurface £ and a collection of 
scalar field to solve the diffeomorphism constraint as it occurs in the second reference 
of I). 

The symplectic reduction of the diffeomorphism constraint has expectedly lead to 
nonlocal 'already smeared' Dirac observables in a natural way. It was not necessary 
to make use of a complete system of mode functions on the hypersurface (exploiting 
the fact that the Hilbert space of square integrable functions on E is separable) to 
write the true degrees of freedom in such a countable fashion. 

The solution of the scalar constraint quadratic in the momenta reveals that the 



29 



Dirac observables will in general mix momentum and configuration variables relative 
to the polarization of the phase space that one started with. 

The formalism of deparametrization applied to field theories shows that there are 
even more serious problems than those that have been dicussed for finite-dimensional 
examples in M at least on the technical side. 

The fact that the loop-variables are so badly suited for our model is quite dis- 
couraging and scratches at the hope that many authors had when employing them 
for gravity, namely that they could serve as the natural 'polarization' to solve quan- 
tum gravity nonperturbatively. Actually, the quantum solutions that have been 
found in |7|] are solutions for the degenerate sector only. In our model, the variables 
T a ( z -) and T*, x \ are simple expressions and the former is easily seen to be a Dirac 
observable on the degenerate sector only. However, it is not the only solution of 
the constraints (in |7j it was speculated that the found solutions are complete). The 
assumption arises that it will be a general feature that the solution of the constraints 
expressed in terms of the loop variables will take a simple form only on the degen- 
erate sector. 

On the other hand, on should never overestimate the results obtained from a model 
and it could in fact also happen that our conclusions are due to an artefact of the 
reduction process. 



A Calculation of the spin-connection 

We will derive here a more general result, namely how the spin-connection looks like 
for the Gowdy models (||). 

The densitized triad for the Gowdy models takes the form 



(3?) 



( Ef EfO 
El Ef 
\ E 2 



so that q 



( EfEf EfE y 
EfE y j E y jE y j 
V (E z ) 2 



(A.l) 



is the expression for the twice densitized inverse metric EfE\. We conclude that 



[det(g)] 2 = [det(^ a )] 2 = (E*det(E^)) 2 =: (E Z E) 2 
and can invert expression (A.l) to obtain the metric tensor 



{(lab) = 



( E y jE y j 



-EfE y j 
-EfE y j EfEf 








V 







o (f) 2 



(A.2) 



(A.3) 



From this we infer the expressions for the triad and its inverse, that is 

E„ a 



e, 



det(g) 



and e l a = 



(A.4) 



namely 



VWe 



( Ef E\ 

E%El I and CO . \\^— 
\ E z 




I E y -El 
-E\ Ef 

V o 







(A.5) 
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which we have to insert into the formula for the spin-connection 



Note that e l a b = {e % a )'5b z in the course of the calculation and that e z e z e l ^ k = e l z {e^)'e l ^ k - 
e\ ef = etc. due to the symmetry properties of the triads so that we can simplify 
(A. 6) without reference to the explicite form of the triads to 



= -\e^[et((eiy + e#(4)') + 5 az e){e k )>\ . (A.7) 



More explicitly 



K = \^{e k x )' + e){e k y )') . (A.8) 

Now we have to go into details and have to make use of the expressions (A.l). The 
calculation is rather tedious. The result is given by 

(K) = ^(^^- + ^[(myE-(mymm + (mymm\, 

+ ±[(E%)'E + {E\)'E X E\ - {El)<E\E%$) 

K) = ^(—^—^ z —^- + -^[—(Ei)'E — (E^'EjEf + (E^)'EjEj}, 
(E Z )'E X 1 

+ e[-(*%ye + {my mm - (my mm], o) 

r z = ^(0, 0, Ej(Ef)' - E^(Ef)') . (A.9) 

These are the corresponding real parts for the Ashtekar connection for the Gowdy 
models. Let us restrict now (A.9) to our model, that is 



E x cos(/3) -Ey sin(/3) \ 
(Ef) = | E x sin(/3) EP cos(P) 

E z J 



(A.10) 



Then (A.9) simplifies tremendously to 



-r x sin(/3) T y cos(/3) 
(r*)= | r x cos(/3) r>in(/3) [ (All) 

o o r 2 

where 

I E y E z , 1 E X E Z , , 

Tx = 2&^~E*~^ Vy = ~2E~^~Ev~^ Tz = ' (A ' 12) 
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